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Localization lengths of superconducting quasiparticles \ s are evaluated and compared with the 
corresponding normal state values \ n in one and two dimensional lattices. The effect of supercon- 
ducting correlation on the localization of quasiparticles is generally stronger in an off-site pairing 
state than in an on-site pairing state. The modification of superconducting correlation to A is strong- 
ly correlated with the density of states (DOS) of superconducting quasiparticles. A s drops within 
the energy gap but is largely enhanced around energies where DOS peaks appear. For a gapless 
pairing state in ID or a d-wave pairing state in 2D, \ s /\ n at the Fermi energy Ef is of order 1 
and determined purely by the value of gap parameter not by the random potential. For the d-wave 
pairing state, the localization effect is largely weakened compared with the corresponding normal 
state and quasiparticles with energies close to Ef are more strongly localized than other low energy 
quasiparticles. 



PACS number: 74.20. Mn, 73.20.Fz 

There has been a renewal of interest in the study of 
the localization effect of quasiparticles in superconduct- 
ing states. Recently in explaining the thermally activat- 
ed behavior of the microwave conductivity in high T c 
superconductors, Lee has proposed that a conductivity 
mobility gap exists in the d x 2_ y 2-w&ve pairing state [1]. 
Since a finite density of states (DOS) exists at the Fermi 
level in a disorder d-wave superconducting state [2] he ar- 
gued that this system should formally resemble a metal 
with a finite Fermi surface where the low energy quasi- 
particles are strongly localized. More recently Hatsugai 
and Lee [3] have numerically studied the effect of impuri- 
ties on Dirac fermions on a square lattice and found that 
low energy states are strongly localized compared with 
high energy states. Their result lends support to Lee's 
argument as the spectrum of the Dirac fermion is similar 
to the d-wave superconductor. However, the supercon- 
ducting correlation is not included in this Dirac fermion 
model. 

The effect of localization in superconductors has actu- 
ally been studied for many years, with the emphases on 
either the change in T c or the effective electron-electron (- 
phonon) interactions [4] . The weak localization of quasi- 
particles in the s-wave superconductor was first studied 
by Ma and Lee [5] in 1985. Within a mean-field ap- 
proximation, they found that superconductivity persists 
below the mobility edge and the quasiparticle excitations 
of the superconductors formed from localized states are 
also localized. Purely off-diagonal disorder (i.e. a spatial 
fluctuated gap function) can also localize superconduct- 
ing quasiparticles. This problem has been addressed by 
Lambert et al [6]. 

In this paper, we study the effect of superconducting 
correlation on the localization of quasiparticles by eval- 
uating the localization lengths of quasiparticles in one 
and two dimensions. The model used here is the weak- 
coupling BCS Hamiltonian with disorder which is defined 
by 



H = t Y. 4aC r >a+J2( A *- C U C l+rl + h - C -) 
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where < rr' > denotes nearest neighbors and fi is the 
chemical potential. V r is a white noise potential which 
distributes uniformly between — V and V. A rr is the 
superconducting order parameter which should be deter- 
mined self-consistently from the Bogoliubov-de Gennes 
equations. Without disorder, A rT is translationally in- 
variant with a particular symmetry with respect to r. 
With disorder the self-consistent gap function A rr is s- 
pace dependent. To solve the Bogoliubov-de Gennes e- 
quations in this case is rather difficult. For simplicity we 
shall ignore the fluctuation of A rr in space. We assume 
A rr = A r and take A r as a free gap parameter. This 
approximation neglects the contribution of off-diagonal 
disorder to the localization of quasiparticles. But as we 
have shown recently [7], the error resulted from the ap- 
proximation is small provided A r is given by the average 
of the self-consistent energy gap A rT . Moreover, in the 
weak disorder limit (V — > 0) the fluctuation of A rr in 
space is negligibly small. Results obtained in this limit 
should be rigorous. 

In the absence of superconducting correlations, H is 
the ordinary disorder Anderson Hamiltonian. It is well 
known that all states are localized in one and two di- 
mensions in this model according to the standard scaling 
theory of localization [8] . The presence of superconduct- 
ing correlation may modify the properties of localization 
of electrons in two aspects. First, the presence of super- 
conducting correlations induces long range correlations 
between electrons which may enhance the coherent mo- 
tion of electrons. Second, the presence of superconduct- 
ing correlations opens a gap in low-lying excitations. The 
opening of this gap suppresses the low energy DOS (it is 
zero within the gap for a clean s-wave superconductor) 
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which may reduce the mobility of electrons and make low 
energy quasiparticles more localized. On the other hand 
the DOS right above the gap edges (E ~ A) diverges 
in the clean limit. With disorder the DOS is no longer 
divergent, but still very large at the gap edge which may 
enhance the mobility of electrons. 

The Hamiltonian (1) is bilinear in fermion operators. 
For the convenience in the discussion below, we perform a 
particle-hole transformation for the down-spin electrons, 
c r i i — > c^, and re-express (1) as 

H' = t ^ (<4t Cr 't ~ C U°r'i) + X]( ArrC rt Cr + T 4. + hx -) 
<rr'> rr 

+ y^(Vr ~ y(i)(<4 t c rt - c^c r4 .) + const. (2) 
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H' has the usual tight binding model form, but the hop- 
ping constant, the chemical potential, and the impurity 
potential have opposite signs for the up spin electrons 
and the down spin electrons. The superconducting pair- 
ing term in (1) now becomes the spin flipping term in 
(2). 

We evaluate the localization length using an iterative 
method applied to very long strips combined with a finite 
size scaling analysis [9]. The localization length on a long 
strip is defined from the one particle Green's function 
G(E) according to the relation 
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where N and M are respectively the length and the width 
of the strip. M=l is the true one-dimensional case. In 
our calculation, N is generally terminated when N is more 
than 100000 in ID or when N is 1000 times larger than 
the localization length A(M, E) in 2D. In averaging A over 
the disorder, we have taken 200 random configurations of 
disorder in ID and 10 random configurations of disorder 
in 2D. The relative error estimated from the deviation of 
A around the mean value is less than 3%. 

We first consider the ID case. An on-site pairing s- 
tate A T =A(5 Tj o and an off-site pairing state A T =A((5 Tj £ + 
&t,-x) a t or away from half filling will be investigated. 
The calculation in ID is simple. It nevertheless reveal- 
s some intrinsic properties of localization effects in su- 
perconducting states, which may help us understand the 
effect of superconducting correlation on the Anderson lo- 
calization in high dimensions. At half filling (fi = 0) the 
off-site pairing state has a gapless excitation spectrum 
with a finite DOS at the Fermi level. This special case 
provides a chance for understanding the localization ef- 
fect in gapless superconductors from ID. 

Fig. 1 shows A s /A„ (the ratio of the localization length 
between a superconducting state and the corresponding 
normal state) as functions of energy E for the on-site 
pairing state with /j, = and the off-site pairing state 
with = and —0.5 in ID. For comparison, the corre- 
sponding DOS is also shown in Fig. 1. Essentially three 



features are revealed by the hgure: 1. A s /A n and DOS 
are highly correlated. For the pairing states with finite 
energy gaps, (b) and (c) in the figure, A s /A„ drops very 
fast within the energy gap. On the other hand, A s /A„ 
shows a peak around an energy where a DOS peak ap- 
pears. The higher the peak in the DOS, the higher the 
peak in \ s /\ n [10]. The drop of \ s /\ n around E ~1.5 in 
Fig lb is a result of the sudden drop in the DOS around 
that energy. 2. For the off-site pairing state A s /A n is 
enlarged in nearly the whole energy region at half filling 
(fi = 0). This enhancement is determined only by A, 
independent of the random potential V, when E is not in 
the region where singular behaviors in the DOS appear 
or very close to Ep (=0 in the figure). For the on-site 
pairing state the enhancement in \ s /\ n is very small out- 
side the energy gap region, very different from the off-site 
pairing state. 3. For the off-site pairing state at half fill- 
ing, \ s /\ n drops near Ep, but only in a narrow energy 
region which is potential dependent (see the inset of Fig. 
la). Low energy states are therefore more localized than 
high energy states in this gapless pairing state. The drop 
of A s /A„ around Ep is not due to the effect of the DOS, 
since in this case the DOS is finite and almost unchanged 
by disorder. Exactly at Ep, A s /A„ is only a function of 
A, independent of V as shown in Fig. 2. This might be 
a common property of gapless superconductors. As will 
be shown later, this universal behavior of A s /A„ at Ep 
appears also in the d-wave pairing state in 2D. 

The difference in the localization effect between the 
on-site and the off-site pairing states is in fact intrinsic. 
To understand this point, let us consider the Hamiltoni- 
an (2). As mentioned previously, after the particle-hole 
transformation, the pairing term becomes the spin flip- 
ping term and the impurity potential has opposite sign 
for up-spin electrons and down-spin electrons. This sign 
difference is not important in the absence of the spin 
flipping term where up-spin electrons and down-spin elec- 
trons are decoupled. However, in the presence of the spin 
flipping term, this sign change has a significant conse- 
quence. Let us consider, for example, a up-spin electron 
hopping from site i, where the potential is negative for 
up-spin electrons, to one of its nearest neighbors, site j, 
where the potential could be either positive or negative 
for up-spin electrons. For the nearest-neighbor pairing 
state, this up-spin electron can hop from site i to site j 
in two ways. One is to hop without flipping spin (i.e. 
hopping through the i-term), the other is to hop and at 
the same time flip its spin (i.e. hopping through the spin 
flipping term, see Fig. 3). If the potential for up-spin 
electrons at site j is negative, the first way of hopping 
is energetically favorable. On the other hand, if the po- 
tential for up-spin electrons at site j is positive, then the 
second way of hopping becomes energetically favorable as 
shown in Fig. 3a. This means that in the off-site pairing 
state the electron can always hop along an energetically 
favorable path by flipping or not flipping spin according 
to the sign of the potential at site j. In other words, it 
means that the effective potential fluctuation < V r 2 > is 
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largely reduced and therefore the potential scattering of 
quasiparticles is largely weakened in the off-site pairing 
state. For the on-site pairing state, however, only the 
first way of hopping is possible. In the case the potential 
for up-spin electrons at site j is positive, this hopping 
costs energy. The second way of hopping can be realized 
through the on-site spin nipping term, but it is not en- 
ergetically favorable as the electron must first jump to a 
higher energy state either at site i (see Fig. 3b, hopping 
process 1) or at site j (Fig. 3b, hopping process 2). This 
shows that the enhancement of superconducting correla- 
tion to the mobility of electrons is weaker in the on-site 
pairing state than in the off-site pairing state in cases 
other conditions are same. 

In 2D, the calculation becomes more complicated. Be- 
low we shall concentrate on a d a .2_ ?/ 2-wave pairing state, 
A r = A(<$ Ti ± x — S Tt ± y ), whose energy gap vanishes at 
some points on the Fermi surface. To find the localization 
length in 2D, the localization length of quasiparticles on 
quasi-lD long strip with finite width M, i.e. \m(E, V), is 
evaluated first. The bulk localization length of quasipar- 
ticles, \(E, V)= liniM hoo \m{E, V) is then determined 

from a single-parameter scaling ansatz 



\m(E,V) 
M 



f A(E,V) 



(4) 



In our calculation, M=7, 11, 15, 19, 23, 27, 31 are used. 
For the d-wave pairing state Am grows extremely fast 
with increasing A or decreasing random potential V. To 
ensure that Am converges fast enough to its infinite long 
strip limit, we limit the study only to small A pairing 
state with relatively large V. 

The scaling function f(x) and X(E, V) are determined 
from a least square fitting procedure [9]. To within nu- 
merical error, we find that the scaling function f(x) (n- 
early log(x) when x is sufficiently large) is the same for 
the d-wave state as for the normal state. This indicates 
that the d-wave pairing state belongs to the same univer- 
sality class, i.e. the orthogonal class, as for the normal 
state, in agreement with the conclusion drawn based on 
the time-reversal symmetry consideration by Balatsky et 
al [11]. On the other hand, it implies that the conven- 
tional scaling theory is applicable and all quasiparticles 
are localized in the d-wave state in 2D. 

Fig. ?? shows \ s /\ n as a function of E for the d-wave 
pairing state. The behavior of \ s /\ n for the d-wave state 
shows a great similarity to the ID gapless pairing state 
in Fig. 1. But the enhancement of A s /A„ when E not 
too close to Ep or to the band edge is much larger than 
its ID counterpart. The enhancement of A s /A„ around 
the band edge is relatively small and nearly potential 
independent, similar to the ID gapless states when they 
are far from the singularities of DOS (the DOS of a clean 
42- s 2-wave state diverges at an energy E ~ A, not at 
the band edge). For the d-wave pairing state, the DOS at 
Ep is largely suppressed by the opening of energy gap. 
But the corresponding \ s is not dramatically reduced 
compared with A„ at Ep. Like its ID counterpart, \ s /\ n 



at (Ep) tor the d-wave pairing state is determined only 
by A, independent of V. It decreases for small A and then 
increases for large A with a minimum located at around 
A ~ 0.25, again similar to its ID counterpart (Fig. 2). 

The large enhancement of A s /A„ in nearly the whole 
energy band in Fig. ?? shows that the potential scat- 
tering of quasiparticles is largely weakened by the su- 
perconducting correlation and the localization effect is 
weaker in the d-wave pairing state than in the normal s- 
tate. The sharp drop of A s /A„ in a narrow region around 
Ep suggests that the states in this region are more local- 
ized than other low energy states. This agrees with the 
results of [1] and [3]. For an infinite layer system with 
each layer a d-wave superconductor, if we allow electrons 
to hop between layers, then those weakly localized states 
may first get delocalized. However, strongly localized s- 
tates around Ep may still be localized when the hopping 
constant between layers t± is sufficiently small. Hence a 
mobility edge may exist at low energy in this quasi-3D 
d-wave pairing state. The value of the mobility gap is A 
and V dependent, but it should be very small as Fig. ?? 
suggested. 

In conclusion, we have studied the effect of supercon- 
ducting correlation on the localization of quasiparticles in 
low dimensions. There is an intrinsic difference between 
an on-site pairing state (for example an isotropic s-wave 
pairing state) and an off-site pairing state (for example 
a d-wave pairing state) as far as the localization effect is 
concerned. The effect of superconducting correlation on 
the localization of quasiparticles is generally stronger in 
an off-site pairing state than in an on-site pairing state. 
A s drops within the energy gap but is largely enhanced 
around an energy where a DOS peak shows. For the gap- 
less pairing state in ID and the d-wave pairing state in 
2D, A s /A„ is enhanced in nearly the whole energy region 
and quasiparticles with energies close to Ep are more 
strongly localized than other low energy quasiparticles 
but not more strongly than the corresponding normal s- 
tates. For these gapless pairing states, A s /A„ at Ep is 
determined only by the value of A. 

I wish thank J. M. Wheatley for many fruitful discus- 
sions. 
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FIG. 1. (a) Ratio of the localization length between a 
superconducting state and the corresponding normal state 
As/A n vs energy E for four white noise potentials for a n- 
earest neighbor pairing state A T = A<$ T= ±i with fi = and 
A = 0.25 in ID. Inset: A s /A s for small E. (b) Same as (a) 
but with fi = —0.5. (c) Same as (a) but for an on-site pairing 
state A,- = A5 T= o. (a'), (b'), and (c') are density of states 
corresponding to (a), (b), and (c), respectively. 




FIG. 2. As/A„ vs energy gap A at Ef for the off-site pairing 
state with fi = in ID. 
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FIG. 3. Schematic representation of the pairing induced 
hopping processes of an electron from site i to site j in (a) 
an off-site pairing state and (b) an on-site pairing state. S- 
cattering potentials are equal in amplitude but opposite in 
sign for up-spin and down-spin electrons at each site after a 
particle- hole transformation (see Hamiltonian H'), which are 
represented respectively by solid and dash horizontal lines. 
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As/An vs E for a d-wave pairing state with /_i = in two 
white noise potentials, V=2 (diamond) and V=3 (cross), 
in 2D. Inset: A s /A n for small E. 
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